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VIBRATIONAL CONTROL OF A NONLINEAR ELASTIC PANEL 


P.L. CHOW* AND L. MAESTRELLO* 


Abstract. The paper is concerned with the stabilization of the nonlinear panel oscillation by an active 
control. The control is actuated by a combination of additive and parametric vibrational forces. A general 
method of vibrational control is presented for stabilizing panel vibration satisfying a nonlinear beam equation. 
To obtain analytical results, a perturbation technique is used in the case of weak nonlinearity. Possible 
application to the other type of problems is briefly discussed. 

Key words, vibrational control, stabilization, nonlinear panel oscillation, perturbation technique 

Subject classification. Applied and Numerical Mathematics 

1. Introduction. The problem under consideration is the stabilization of the nonlinear panel oscillation 
by an active control with a vibrational actuator. This work was motivated by the recent experimental 
investigations of the second author (Maestrello [1, 2]), who demonstrated clearly that the vibrational control 
could be an effective means of stabilizing the boundary-layer flow as well as the panel vibration. This paper 
will offer a general method of vibrational control and its application to the problem involving a nonlinear 
elastic panel excited by the periodic wall-pressure fluctuation in a boundary-layer flow. 

The general principle of active control, the vibrational control in particular, is to introduce an action 
which affects a change in the behavior of a dynamical system in a desirable manner. In the boundary layer 
transition control [1], the periodic heating and cooling of the wall induce a parametric vibration of the fluid 
viscosity which, in turn, stabilizes the flow. In the case of panel vibration [2], a properly added vibrational 
force with the same forcing frequency may result in suppressing the subharmonic oscillations (see Section 
2). The suppression of subharmonics, of course, has the implication of controlling the chaotic motion. 

The main idea of vibrational control stems from the fact that an inverted pendulum can be stabilized 
at its upper equilibrium position when the lower suspension point executes a rapid vertical vibration (see 
e.g. [3]). Based on this idea, a general principle of vibrational control was proposed, notably by Meerkov [4], 
to stabilize the equilibrium points of some finite- dimensional linear systems. Application of this principle 
to reactor dynamics was done by Bellman et al. [5]. By contrast, in this paper, we will extend this control 
principle to stabilize the periodic motions of infinite- dimensional systems, instead of equilibrium points in 
finite dimensions. In addition to the high-frequency parametric vibrational control used in [4, 5], a vibrational 
force with the same forcing frequency will be required. Unlike the usual feedback of feed-forward control, 
the vibrational control does not need accurate measurement of the system inputs and outputs and can be 
implemented much more easily, especially for an infinite- dimensional system under consideration. 

In this paper we consider the panel vibration which satisfies the initial boundary value problem for the 
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nonlinear beam equation [6]: 


( 1 . 1 ) 


mdfw + cd t w ~{Q + N(t)]d*w + Dd*w = A p(t, x), 0 < x < l 
< w(t , 0) = w(t , l) = 0, <9 2 tu(t, 0) = d%w(t, l) ~ 0, 

, w(0,x) = w 0 (x) y d t w( 0,x) = (x). 


Here w denotes the transverse deflection; d t ,d x ... are the partial differentiations in f,x...; the positive 
constants m,c and D represent the unit mass, the damping coefficient and the bending stiffness of the panel, 
respectively. The axial force Q is positive or negative according to the force being tensile or compressive. 
The large panel deflection introduces an additional tension N(t) given by 

(1-2) N(t) = bf |<9 x u;(t,x)| 2 dx, 

Jo 

where b is an elastic constant. The forcing term A p denotes the pressure difference across the panel surfaces. 
The homogeneous boundary conditions mean that the panel is simply supported, and the initial data w 0 
and w\ are given. Suppose that, without any control, the periodic solution of equation (1.1) excited by the 
pressure A p is unstable. Our problem is to stabilize the panel oscillation by applying an appropriate control 
in the form of vibrational forces added to the axial force and the pressure A p. 

The paper is organized as follows. To illustrate the basic ideas involved, in Section 2, we consider the 
control of the Duffing equation, for which the response characteristics to a time-harmonic excitation is well 
known. The feasibility of the vibrational control can be discussed geometrically by referring to the response 
curves. Since the applicability of vibrational control is not limited to the structure dynamics, in Section 
3, a general method of vibrational control for a class of nonlinear evolution equations is presented. For a 
given unstable periodic solution, the control strategy is to shift the Liapunov exponent r of the vibrational 
equation to the negative half-line so that the corresponding periodic solution becomes stable. This method 
is applied to the nonlinear panel vibration problem satisfying equation (1.1). For weak nonlinearity, analytic 
results are obtained by a perturbation analysis and the case of single-mode excitation is worked out in detail. 
Finally, in Section 5, some concluding remarks are made and other possible applications such as the flow 
stabilization problem are mentioned. 


2. Control of Dufflng’s Equation. Before dealing with the nonlinear beam equation (1.1), we con- 
sider the Duffing equation 


( 2 . 1 ) 


y + fiy -f Sy + fiy 3 = F cos u;t, 


where the dot denotes the time derivative, the constants jU, and F are ass um ed to be positive here, and 
w > 0 is the forcing frequency. For small F, by perturbation analysis [7], it is known that equation (2.1) has 
a periodic solution of the form 


( 2 . 2 ) 


y — A cos (u>t + 0) 


for some phase shift 0, where the amplitude A is related to the frequency t o by the response equation [8]. 


(2.3) 


{uj 1 - 5) A - -BA 3 ] + iiurA 2 = F 2 . 
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Fig. 2.1. Response curves for harmonic oscillation 

By varying the value of F, equation (2.3) yields a family of response curves in the \A\ — oj plane. Referring 
to Fig. 2.1 for F = Fo, Fi, the solid portion of the curve corresponds to the stable regime for the periodic 
solution (2.2), while the dotted part of the curve (between two points of vertical tangency) renders the 
solution unstable. With the frequency uj fixed, point U on the F 0 -curve is unstable, but, by changing Fq to 
Fi, point U moves up to point S on the Fi-curve becoming a stable point. Therefore, in this case by adding 
an in-phase force with the same frequency, an unstable periodic motion can be stabilized. On the other 
hand, if the forcing amplitude F is large, the system may exhibit a subharmonic response. For example, 
consider the case of subharmonic response with frequency Again, by a perturbation analysis, it is found 
that equation (2.1) has a subharmonic solution of the form [8]: 

y = A cos (cot + 6 i) 4* B cos -f 62 j 

where A^B and 0 i ,02 are the corresponding amplitudes and phases, which satisfy some response equations. 
For the subharmonics, the equation reads 

(T)-f 

with / = 9F/8. For F = F 0 ,Fi with Fi > Fo > 0, the response curves associated with (2.4) are given in 
Fig. 2.2. Note that for u; = u?o, point P on the Fo-curve corresponds to a subharmonics with amplitude 
ro. However, this subharmonics will disappear when F changes from Fo to Fi, since lj < u>o, and u)\ is the 
smallest frequency for the existence of a subharmonics at F — F\. This may explain qualitatively why a 
subharmonic vibration can be suppressed in the experimental investigation [2] by an additive periodic force, 
which has the effect of changing the forcing amplitude F. 

In contrast with the additive vibrational control, the control can be applied parametrically. For instance, 
we regard the Duffing equation (2.1) as an approximate equation for an inverted pendulum near the upper 
equilibrium position ( y = 0), for which a = — S < 0. Clearly y = 0 is an unstable equilibrium. If the 
suspension point vibrates at a high frequency v the equation (2.1) should be replaced by [3] 

(2.5) z -h fiz + \p(vt) — S] z + flz 3 = F cos ut, 


(2.4) 


U ) 2 = 9a + ^7 ( B 2 + 2 f 2 ) ± 
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Fig. 2.2. Response curves for subharmonic oscillation 


where p(r) = p(r -f 2tt) is a periodic function. Without the control p, a periodic motion about y = 0 is 
obviously unstable. However, by the method of averaging [3], equation (2.5) can be closely approximated by 
the averaged equation (see Section 4) : 


where < p >= 0 and 


y + py + [< p 2 > - y + /3y 3 — F cos cut , 


1 f 2 n 

< p n >= — / p n (r)dr, for n = 1, 2. 
27T J 0 


Therefore, if a =< p 2 > — S > 0, the periodic motion can now be stabilized as before. 

The above examples show the possibility of stabilizing periodic motions by vibrational control. As a 
generalization we consider the following control problem: 

(2.7) z 4- pz -f az + flz 3 = /(A, a ;t) + h(A, vt , z) 

where /(A, r) = /(A, r + 2n) and h( A, r, z) = /i(A, r + 27r, z) are periodic functions; A is a control parameter 
and h is a certain control function with /i(0,j4,y) = 0 and uj>v are the vibration frequencies with v oj. 
The uncontrolled case corresponds to A = 0 and h = 0. Of course the equations (2.1) and (2.5) are special 
cases of (2.7). Suppose that z == is an unstable periodic solution of equation (2.7) when A = 0 and 

h = 0. The control objective is to choose the control parameter A and Function h so that the corresponding 
periodic solution z = <p(A, t) with <p(0 ,<) = <po(£)> becomes asymptotically stable. This means analytically 
that the variational equation for y = (z — <p) from (2.7) has only exponentially decaying solutions. More 
precisely, if 

r(X,ip) = lim jin |j/(t)|, 

t—>oc t 

then we must choose A and h such that r < 0. Obviously, unlike the optimal control, such a control, if 
possible, is far from unique. The choice of A and h, though guided by physical feasibility, is mostly up to the 
personal preference. In what follows, this control principle will be generalized to deal with nonlinear partial 
differential equations. 
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3. Stabilization of Nonlinear Evolution Equations. In the theoretical discussion, it is convenient 
to consider the partial differential equations of interest as a nonlinear evolution equation of the form 


(3.1) 


= B{u) + F( A, ut) + H( A, ut , u), 

u(0) = 


where n(t) is a vector in some infinite-dimensional vector space V with initial state h. The operator B is 
nonlinear, F(A, r) = F(A, r 4- 2n) and H(X i r ) u) — H( A, r -j- 27r, u) are periodic with control parameter A, 
being a scalar or vector. The control function H acts parametrically with rapid oscillations so that v >> uj. 
When A = 0, H(0, r, u) = 0 and the system (3.1) is uncontrolled. We are interested in stabilizing a periodic 
motion which is unstable at A = 0. If the equilibrium solution uo of (3.1) at A = 0 is also unstable, we 
introduce a parametric vibrational control H to stabilize it as in the case of an inverted pendulum. The 
effect of H can be examined by the method of averaging [3]. By a change of time from t to a = tjs with 
£ = l/v, the system (3.1) can be approximated by the averaged equation 

f % = B(u) + F(\M) + H(\,u), 

\ «( 0 ) = h, 

where 


1 f 2n 

(3.3) H( X, u) = (H(X, ut,u)) = — H( A, r, u)dr. 

Jo 

The function H should be chosen so that the equilibrium solution of Ui of the averaged equation (3.2) 
becomes stable. Without control (A = 0 ), let u = be an unstable periodic solution of equation ( 3 . 1 ) 

near ito* In addition to the parametric control H , we have modulated the forcing function F(A, u>i) by tuning 
the control parameter A so that the corresponding periodic motion satisfying the averaged equation ( 3 . 2 ) is 
asymptotically stable. To this end let us consider the variational equation of ( 3 . 2 ) for v = (u — tp): 


(3.4) 


where 


dv 

dt 


B i(ip,v) + Hi(X,ip,v), 
v{0) = 9, 


Bi(ip,v) = B(v + ip) - B(ip), 
Hi{X,ip,v) = H(X,v + i>)-H( X,rP), 


and g is an initial vector in V. Let \\h\\ denote the magnitude (norm) of vector h. The control objective is 
then to choose function H and parameter A in such a way that the Liapunov exponent r is negative, 

(3.5) r(A, H) = lim hn ||n(t)|| < 0 , 

t~* OO t 

for all g with j [ < 7 1 1 < 5 with some S > 0 . For small <5, the variational equation (3.4) can be linearized to give 


(3.6) 


r %=A(\, t )v, 

\ u(0) = 9, 


where A(X,t) is a linear operator defined by 


(3.7) 


A(A, t)v = {B u [V>(A, t)] + H u [A, V’( A, t)]} v, 
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and B u (v) = S - B ^ ,H U { A,t>) = are the linearized operators of B and H at v. Note that A(X,t) is 

periodic with the same period T as that of ip. *>,w, if V is finite-dimensional and A( A, t) is a matrix, then, 
by the Floquet theory [7], the solution of equation (3.6) can be expressed as 

(3.8) v(t) = P{t)e iR g, 


where P(t) = P(t -f T) is a periodic matrix, and R is a constant matrix. The smallest real part of the 
eigenvalues of R yields the Liapunov exponent r. Of course the representation (3.4) holds for any finite- 
dimensional approximation of equation (3.4). Unfortunately, even the periodic function ip is known, the 
analytical computation of the Liapunov exponent r through either (3.5) or (3.8) is impossible without 
simplifying assumptions. For example, for small amplitude vibration, the nonlinearity is weak so that the 
perturbation method and an eigenfunction expansion can be applied. This procedure will be illustrated in 
the application to the panel vibration problem. 

4. Vibrational Control of Elastic Panel. By redefining the constants in the nonlinear beam equation 

( 1 . 1 ) under a vibrational control, it yields 

(4-1) d 2 t w + (j,dtw — (a 4* /3\\d x w\\ 2 ) d x w + 7 d*w = p(\,u?t, x) + h( A, ut, x, w) 


where the initial-boundary conditions are omitted, and 


(4.2) 



(4.3) p( A, u>t , x) = Ap(ut, x) + pi (A, u )t, x), 

pi and h are the additive and parametric control forces with frequencies u u. The physical constants 
p, /3 , 7 arc positive, while a is positive or negative depending on the axial force being tensile or compressive. 
Without control, we assume that, at A = 0 ,pi( 0 , x, r) = h(0,a, x, w) = 0 . To be specific, we choose the 
parametric control to be a vibrational axial force of the form 

(4.4) h{ A, ut , x, w) = g(A, i/t)d%w, 
where 

q = d t q{\, ut) 

with 

(4.5) < q >=< q >= 0 
Now let u\ — w and u 2 defined by 

(4.6) ii\ = ^2 + q( A, ut)d 2 u\ 

Then equations (4.1) and (4.6) yield 

ui = ~ {v( u 2 + qdlu x ) - (a + (3\\d x ui\\ 2 )dlux + 

V " ' 2 - qd£u 1 ) + p(X,uit,x). 


;i 1 r 
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We set 


u = 


Ui 

U 2 


and rewrite the equations (4.6) and (4.7) in the form (3.1): 


(4.8) 
where 

(4.9) 

(4.10) 
and 

(4.11) 


B(u) = 


du 

— — B(u) + F( A, u it) + H( A, ut , tx), 


U2 

{nm - (a + P\\d x u i|| 2 ) dlui +7^ui} 


F( A, cjt) 




i/(A, i4, u) = 


Qdx u i 

-fiqd^ui - q 2 d*m - qd 2 u 2 


In view of equations (4.5) and (4.9) (4.11), by taking the time-average of equation (4.8) in a = ut with 
r = ut fixed, we get 


(4.12) 
where 

(4.13) 


du '■ 

— = B{u) + F(A, ut) + H( A, tx), 


H{X,u) = 


0 


- <q 2 > d*u x 

We note that the average equation (4.12) yields a scalar equation for w = Ui as follows, 

( 4 - 14 ) dfw + fidfw - (a + 0\\d x w\\ 2 ) d\w + ( 7 + < q 2 >)w = p(X,ut,x), 

which shows that the high-frequency axial vibrational force q(X ) ut) has the effect of increasing the bending 
stiffness 7 by the magnitude of < q 2 >. Thus it stabilizes the system statically in general. Now let 


^oi(£) 

*pQ2(t) 


uo = ^0 (0 = 

be an unstable periodic solution of equation (4.8) when A = 0 (H - 0), and let 

u — ip(X, t) — 


i>i (x,t) 

^2 (A ,t) 


be a periodic solution of equation (4.14) with ip(0,t) = ^o(0- Define v = (u - ip) so that v satisfies the 
variational equation (3.4). Here it can be written in the form 


(4.15) 


dv 

dt 


— A(X, t)v -j- j3G{ A, t, u), 
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where A is a periodic linear operator and C, is a nonlinear mapping defined as 


Av = " , G = 

h 


fi = ~nv 2 + (a+ < q 2 > +/?||d x y>i|| 2 ) d%vi +2P(d x ipi,d x v 1 )d^ ] . - 7 ^ 1 , 
f 2 = ||9x^i|| 2 ^^i + 8^)81, + \\d x vi\\ 2 3%vi, 


with the inner product notation 


(g, h)= g(x)h(x)dx. 
Jo 


When the nonlinear term G is dropped, equation (4.15) yields a generalized Hill’s equation, a linear partial 
differential equation with periodic coefficient: 

(4.16) ^=A{\,t)v. 

For computational purposes, introduce a complete set of orthonormal functions {e n }, which may be the 
eigenfunctions associated with the linearized problem, or e n (x) = ^/f sin ^f-x^n — 1 , 2 , ... By the expansion 
of the solutions of (4.12) and (4.15) into terms of e ; n s as follows, 

oo oc 

u = 52 u n (t)e n ; v = 52 v n(t)e n , 
n= 1 n= 1 

their coefficients satisfy the infinite systems of coupled ordinary differential equations of the form: 

du ■ * 

i it = Bi(m , ..., Ufu . ..) H" F t (A, £) “h Hi (A, Wi, u n , •••), 
at 


UVi V”^ / \ \ . 

^ (A, 5 ^ 1, 2, ... , .. ., 

n=l 

where a t j = (Aei,ej). The above systems can only be solved numerically for truncated systems of low 
dimensions. As mentioned before, if the nonlinear effect is weak, we can apply the perturbation analysis to 
approximate solutions analytically. To this end let us assume that the damping coefficient and the forcing 
amplitude are small. By proper scaling with respect to a small parameter e > 0, the equation (4.14) is 
rewritten as 


(4.17) 


dfw + efidtw — (a + ef3\\d x w \\ 2 ) d\w -f 7 d\w = £p(A,wt,x), 


for which we assume a > 0 and set < q 2 >= 0 for simplicity. It remains to study the problem of additive 
control. To illustrate the perturbation procedure, we will analyze the case of single-mode excitation in some 


detail. 


Let us consider the case of nth mode harmonic excitation in (4.17): 


(4.18) 


p(A,utf,x) = F(A) sin — xcos ujt^n = 1 , 2 ... 


where the control parameter A modulates the forcing amplitude F and F 0 — F( 0) is the uncontrolled 
amplitude. Then the equation (4.17) admits a single-mode solution 


(4.19) 


. , nn 

w — z n (t) sin —x 

l 
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Fig. 4.1. Response curves for different modes of excitations 


and z n satisfies the Duffing equation: 

(4.20) Zn + + ct„z n + = eF( A) cos cut, 

where 

a„ = [a + 7 m 2 ]m 2 . 

(3 n = f(^) 4 - 

The perturbation analysis of Duffing’s equation has been discussed by many authors (see e.g. [3, 7]). Here 
we adopt the method of averaging by letting 

(4.21) Zn = Vi (t) sin u)t 4- y 2 (t) cos uit 
with 

yi sin wt 4- 7/2 cos cut = 0 


which are substituted into (4.20) to give 

(4.22) 


V\ = 2 b [ S nV2 - \(3n\y\ 2 V2 + + ^O)] , 

i/2 = ^ [<Wi - \Pn\y?y\ +y-wy 2 ] , 


with S„ = (w 2 — a n ) and \y\ 2 =y\ + y 2 - In the polar form, y\ = r sin if, y 2 = r cos if , this equation becomes 

(4.23) 
where 


r = <P,u), 

<P= 


(4.24) 


R — fiwr — F sin ip 

$ = 6n- \Pn r 2 + f COS If. 


Z n = T(t) COs[a;t + 


The solution (4.21) can be written as 
(4.25) 
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Therefore, for z n being periodic with frequency ui,r and tp must be constants, which correspond to the 
equilibrium point (f,tp) of equation (4.23) satisfying 


(4.26) 


R(r, tp, uj) = 0 

<f>(f, <£, w) = 0. 


By taking equation (4.24) into account, the above equation can be solved approximately to give 
(4-27) 2 + M f 2 u, 2 = F 2 ( A), 

which, by a change of notation, agrees with the response relation (2.3). Therefore for each n, the response 
curves are shown in Fig. 2.1. Schematically, for n = 1 , 2 ,..., the response curves are plotted in Fig. 4.1. 
Geometrically the control strategy is to steer an unstable point I7 n on the F 0 -curve to a stable point S n on 
the Fi -curve. Analytically the stability of a periodic solution is now reduced to that of an equilibrium point, 
which can be checked more easily. To do so we form the first variational equation of (4.23) about (r, <p) — 

(4.28) 

where Rr = d r R(r,ip,u>), R^ = d^R(r,(p,u)) and so on. Let 77(A) denote an eigenvalue of the coefficient 
matrix of (4.28). It can be readily verified that, by making use of (4.24) and (4.27), if 

( 4 -29) D{ A) = - f^ r ) > 0, 

then Re 77(A) < 0 so that the steady state (f, <p) is stable. This is of course the stability condition for the 
associated periodic solution. After computing the partial derivatives in (4.29), it yields 

(4-30) D( A) = ^/3 n f 2 - 6 ^ - J n ) + /x 2 u/ 2 > 0, 

with 5 n = (a; 2 — a n ). The above inequality determines the stability regime S in the f - u) plane. In view of 
(4.27), r(A) depends on the control parameter A, which will take an unstable point into the stable regime S. 
Note that from (4.30), we can get a simple sufficient stability condition: 


P 2us f^rP d” i 

® = 2^ (®r/> + R^r^j , 


(4.31) 


3 

u> 2 < a n + -0 n r 2 A, 

4 


or 


w 2 > a n + ^/? n f 2 A. 

4 

The above inequalities give rise to stable (shaded) sub-regions as shown in Fig. 4.2. 

In general, all modes are excited by periodic pressure fluctuations. For instance, consider the harmonic 
forcing (4.18) with a general spatially dependent amplitude 

(4-32) p(X 1 u;t i x) = F(A, x) cosu;t. 


If the axial load is compressive (a < 0) and slightly exceeds the lowest buckling load 7(f) 4 , the parametric 
control Q{yt ) can still be used to stabilize the system statically by choosing <T <7 2 >> |cif|. So we remain to 
consider equation (4.17) by assuming a > 0 there. To apply the above perturbation procedure, we need to 
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nr t r; 



Fig. 4.2. Stability regions (shaded) for different modes of excitations 


expand w in (4.17) and F in (4.32) into infinite series with respect to the modal function e n (x) 
for n = 1,2,... The resulting infinite system of coupled nonlinear differential equations for 
functions can then be treated by a perturbation analysis. Such a procedure developed previously for nonlinear 
wave equations by one of us (Chow [9]) can be applied here. However, unlike the single mode situation, 
simple stability conditions such as (4.30) or (4.31) are no longer attainable. Though it is possible to study 
the stability regime numerically after a finite- mode approximation, this has not yet been done. 



the coefficient 


5. Concluding Remarks. In the paper we present a general method of vibrational control for a 
certain class of nonlinear evolution equations with a particular reference to the nonlinear beam equation 
arising from the panel structure dynamics. The control consists of a high frequency parametric vibration 
and the forcing amplitude modulation. The high-frequency control is to affect a change in system parameter 
for static stability, while the additive control of the excitation force, if needed, is to stabilize an unstable 
periodic motion. In application to the panel structure, we show that, for a periodically excited panel near a 
buckled state, a high frequency oscillatory axial force can keep the system in the state of periodic motion, 
which can then be stabilized by an additive force modulation. The reason that we only control the force 
amplitude, instead of both the amplitude and phase is that the additive control is the most effective when it 
is in phase or out of phase with the excitation force. For a small forcing amplitude, a perturbation technique 
can be used to reduce the stabilization of a periodic motion to that of an eq uili brium point, the latter of 
which is much simpler to analyze. In the case of a single-modal approximation, the stabilization problem 
can be studied numerically but has not yet been treated. The vibrational control principle described in this 
paper can also be applied to other problems such as the flow stability control. Here the nonlinear evolution 
equation is given by the Navier-Stokes equation. For a slightly unstable flow, the perturbation analysis by 
Keller and Kogelman [10] can be employed to deal with the flow stabilization by vibrational control. 
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